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Abstract 

To model ferromagnetic material in finite element 
analysis a correct description of the constitutive re- 
lationship (BH-law) must be found from measured 
data. This article proposes to use the energy den- 
sity function as a centrepiece. Using this function, 
which turns out to be a convex function of the 
flux density, guarantees energy conservative mod- 
elling. The magnetic field strength can be seen as 
a derivative with respect to the fiux density. Es- 
pecially for anisotropic materials (from lamination 
and/or grain orientation) this method has advan- 
tages. Strictly speaking this method is only valid 
for anhysteretic and thermodynamically stable ma- 
terial. 

1 Introduction 

This essay discusses one aspect of numeric mod- 
elling of magnetic fields: anisotropic magnetic ma- 
terial properties. In practical electromagnetic appli- 
cations anisotropic materials are common as lam- 
inated electric steel with and without grain orien- 
tation. In the design process finite element analy- 
sis (FEA) is widely used and correct modelling of 
the material properties is essential. The literature is 
wide spread, but considering saturation effects this 
essay contributes modelling that does not conflict 
with energy conservation. 

Application with use of electric steels are trans- 
formers, motors, and generators. The main reason 
for lamination is that eddy currents in stacking di- 
rection are suppressed. This lamination also has the 
effect that the magnetic properties are anisotropic. 
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Two different types of electric steel are used: grain- 
oriented and non-grain-oriented steels. The non- 
grain-oriented sheets have magnetic properties that 
are isotropic in the sheet plane, whereas the non- 
grain-oriented steels have different magnetic prop- 
erties in rolling direction and in transverse direc- 
tion. 

The discussion here is only applicable for perfect 
soft material, i.e. material which has a nonlinear 
BH-characteristic but has no hysteresis (meaning 
that the area of the hysteresis loop is empty). This 
is the requirement that there exists a one-to-one 
relationship between the vectorial flux density B 
and the vectorial magnetic field H. Models for hys- 
teresis (Jiles-Atherton, Preisach, and other) are not 
considered here. 

The centrepiece of this essay is the use of the 
magnetic energy as a function of the flux density as 
the centrepiece of the description of magnetic prop- 
erties. The magnetic field strength is in this concept 
the gradient (with respect to the flux density) of the 
energy density. The energy density can be used to 
describe linear and nonlinear as well as isotropic 
and anisotropic material. It will be used to derive 
equations for lamination effects and as a base for 
the interpolation rule that extends measured data 
of grain-oriented steel to a full 3-D description. 

Also the energy density function will be used to 
graph BH-characteristics by plotting lines of equal 
energy (iso-lines). 



2 Maxwell-Equations 

In order to clarify the notation, the relevant equa- 
tions for magnetostatics are repeated. The flux den- 
sity is described by a 3- vector B which is divergence 



free (Gau6-law): 



V-B = 0. 



The material properties are in the form of a func- 
tion that maps B to the magnetic field H which is 
also a 3-vector. Formally: 



H = F(B) . 



(1) 



Ampere's -law relates the magnetic field to the elec- 
tric current density vector j, such that the curl of 
the magnetic field is equal to the current density: 

VxH=j. 

Many of the commonly used FEA software uses the 
vector potential A as the unknown function, which 
is linked to the flux density by 

B = V X A 
= V-A. 

The topic in this text is to discuss anisotropic ma- 
terials which are described with a general nonlinear 
function F. 

The differential reluctivity v is a 3-by-3 matrix 
and defined as the partial derivative of the field 
strength with respect to the fiux density. The com- 
ponents are 

_ dHj 

3 Literature review 

Other authors often express the constitutive rela- 
tion by employing the effective reluctivity matrix 
(or equivalently the effective permeability) in the 
form H = i/e//B, where Veff in general depends 
on the field quantities. 

The case of lamination is often treated in the fol- 
lowing way (eg. [H]): assuming steel has the mate- 
rial law 
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and the lamination fill factors are /i for insulation 
and /2 for steel. The BH-relationship for the com- 
posite is 
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But this relationship is only correct if the reluctiv- 
ity is constant. In reality the reluctivity must de- 
pend on the flux density and the reluctivity func- 
tion of the underlying steel must be evaluated in 
the model using the flux density in the magnetic 
material (and not the macroscopic value). The ref- 
erence [7] gives a correct model, where one addi- 
tional equation must be solved. This agrees with 
the solution presented here, but is only valid for 
isotropic material. 

When modelling the in-plane anisotropy of grain- 
oriented electric steel often only a limited num- 
ber of sets of measure data (usually in rolling and 
transverse direction) is available. Then interpola- 
tion models must be found for intermediate direc- 
tions. Usually a 2-D model is often presented in 
literature with the constitutive relation 
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Such an effective reluctivity is not uniquely de- 
fined. There are two physically measurable degrees 
of freedom (ff^j Hy )^ but the symmetric reluctivity 
tensor contains three values {vxxii^xy = i^yxT^yy)- 
So when determining the reluctivity tensor a gauge 
has to be chosen, but no author actually discusses 
the implication of such a choice (eg. [3l fTOl [2l |4]). 
Without going into the detailed equations the 
models all represent BH-relationships with non- 
symmetric differential reluctivity, and therefore vi- 
olate energy conservation as detailed below. 

In cases when measurements are available in the 
entire B^;, i?j,-plane the interpolation can be done 
directly on this data (0). 

The energy density function (or co-energy den- 
sity function) is only used by few authors (eg. 
[T^ [S])- The results presented here can be seen 
as en extension of these ideas so that the en- 
ergy density forms the centrepiece of all treatment 
of BH-relationships and with special attention to 
anistropy from lamination and grain-orientation. 



4 General ansiotropic BH- 
functions 

In the introduction the BH-function (eq. [T]) was 
considered to be a general vector function. Now the 
properties of this function are discussed. 
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Let 7 be a path in the space of B, i.e. a smooth 
function [0, 1] — >■ i?^. The energy density w related 
with this path is defined as 

w{-f) = f H • j'dt. 
Jt=o 

This energy is stored in the magnetic field at one 
material point. Energy conservation now means 
that w depends on the start and end points alone 
and not on the path. Using zero as the start point 
the energy density can be defined as a function of 
B without referring to a particular path: 

w{B) ^ [ H-dB. 

Jb'=o 

The following statements follow from the existence 
of such an energy density function. 

• The magnetic field is the gradient of the energy 
density 

H = VbW 

• Since the differential reluctivity tensor is the 
second derivative of w and is symmetric 

_ dHi _ d^w 
" ^ dB.dBj 

must be fulfilled. 

5 Convex energy function 

For isotropic BH-characteristics it obvious that |H| 
is monotonously rising with |B|. In consequence the 
derivative (the differential reluctivity) is positive. 
Ref. [8] show that for a thermodynamic stable ma- 
terial the differential permeabilty (and its inverse 
the differential reluctivity) must be positive. In an 
anisotropic case the reluctivity is a symmetric ma- 
trix and the requirement for positiveness is gener- 
alised as follows. A matrix M is called positive def- 
inite if for all vectors v the product v ■ Mv is always 
positive. It is therefore required that the differential 
reluctivity is such a positive definite matrix. 

The requirement for a positive definite second 
derivative can be linked with the term of a convex 
function. A scalar function / : i?" — >■ i? is called 



Magnetic energy density for vacuum 




Figure 1: Contour lines of the energy density func- 
tion for vacuum and direction of magnetic field 
strength 

convex if for all x, y e i?" and all t G [0, 1] the 
inequality 

f{{l-t)x + ty) > {l-t)f{x)+tfiy) 

holds. In words: the straight line connecting two 
points in the graph of the function lies above the 
function. The word convex can be explained, when 
looking at point sets like 

U{Q) ^{xe R"\f{x) < Q} (2) 

for Q e f{R"). For convex functions these sets are 
convex in the ordinary sense. A set U{Q) includes 
all points where the function takes values smaller 
or equal than Q. An appropriate way of displaying 
such functions is to look at iso-lines (where f{x) = 
Q), because theses iso-lines enclose convex areas or 
volumes. 

It can be demonstrated that for moderate ad- 
ditional assumption that convex functions have 
positive definite second derivatives and vice versa 
(see n). 

The requirement for a positive definite reluctivity 
is also mentioned in [T3]. It has the benefit that it 
guarantees a stable finite element analysis. 

6 Examples 

For simple cases the energy density functions are 
as follows. A natural way of graphing the energy 
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Figure 2: Contour lines of the energy density func- 
tion for a linear anisotropic material. 



Figure 3: Contour lines of the energy density func- 
tion for mild steel. 



function in 2-D planes is to use contour lines of 
equal energy. These lines enclose convex areas if 
the energy function is valid. 



In vacuum: 



B2 

2^0 



(3) 



with /io = 47rl0^^ Vs/Am. As shown in fig. [l] 
lines of constant energy density are concentric 
circles around zero. This energy leads to the 
magnetic field strength: 



H 



B 

Mo ' 



(4) 



The mentioned figure also displays the direc- 
tion of the field strength H and shows that the 
field strength is perpendicular to isolines of the 
energy density. 

Anisotropic linear case: 

w(B) = B • i^B, 

where ly is the symmetric reluctivity matrix. 
The contour lines are concentric ellipses (fig. [2]) 
and the BH-law is 

H = lyB. 



Nonlinear isotropic case: the energy density 
only depends on the modulus of the flux den- 
sity 



w(B) 



|B|). 



(5) 



The contour plot is similar to the plot for vac- 
uum, only the spacing between the lines is dif- 
ferent (fig. [3]). The underlying BH- relationship 
can be measured and is displayed in fig. [s] The 
resulting BH-rclationship is 



B 



(6) 



The derivative 



is exactly the result of 



measurements of BH-curves. 



7 Lamination energy function 

In principal, modelling a laminated core would re- 
quire to resolve the individual layers of magnetic 
material and non-magnetic insulation in the FEA, 
which is due to the high number of sheets imprac- 
tical. Therefore the method for homogenisation is 
used which models effective B and H fields instead. 
To derive the model the following is assumed: 
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• The coordinate system is such that the sheets 
are in xy-plane and stacked in z-direction 

• Material 1 has the properties of vacuum (cq. [S]) 

• The magnetic properties for material 2 is de- 
scribed by a general energy density W2{'B). 

• The thickness of the sheets is so small that 
changes in the field from one sheet to its neigh- 
bouring sheet can be neglected, and in the 
same way from one insulation layer to the next 
(the fields in the sheet and the insulation, how- 
ever, can be different) 

The volume ratios material 1 and 2 are called /i 
and /2 respectively. 

From basic Maxwell-equations it is known that at 
interfaces the normal components of the flux den- 
sity is continuous (Quantities with indexes 1 or 2 
refer the respective materials, quantities without 
this additional index are homogenised values): 



Biz 



Bo 



Homogenised flux density components are intro- 
duced for the non-continuous components by us- 
ing a mixing rule, which is the natural choice from 
Maxwell-equations 

Bx '■— fiBix + 

By ■= flBiy + /2-B2y 

Now, the homogenised energy density is a mix- 
ing of the energy density of the sublayers. Using 
Bix,Biy as the yet unknown flux densities in the 
insulation layer the energy density can be written 

as 

Wlam [Bx ,By,Bz, Bix , Biy) = 

4^SBl^ + B% + Bl) 
+/2W2 (^j^i^x — fiBix), j^iBy — fiBiy),B. 

which is the weighted sum of the energy densities 
in the two types of material. The flux density vec- 
tor components are the values inside the respec- 
tive material. Minimising the energy with respect 
to Bix, Biy gives 



w/am(B) = min Wi^^{Bx, By, B^, Bix, B 



ly) 



determines the unknowns Bix,Biy and fixes the 
energy density. Evaluating the partial derivatives 
and equation them to zero, results in the equations 
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- Fx{j^{Bx - fiBix), j^{By - fiBiy),Bz) 
= Fy{-h{Bx - fiBix), j:{By - fiBiy),Bz), 



which constitute a system of nonlinear equations 
for Bix, Biy. Deriving with respect to Bx, By, B^ 
defines the magnetic field H. 
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Putting the equations together the following set 
of equation describes a stack of the laminated elec- 
tric steel: 



Hx 
H,, 



/2 



fiB, 
Mo 



{Bx — fiHxf^o) 

= F I j^iBy - flHyflo) 

B, 



Given the effective flux density the magnetic field 
can be computed by solving this system of non- 
linear equations. This is now the material law for a 
stack of electric steel, if the underlying material law 
of the steel type is known (which will be discussed 
below). For general functions F these equation can- 
not be solved analytically. The system of equation 
has to be incorporated into a FEA software in a 
way that, every time the system enquires the value 
of the magnetic field for a given flux density, the 
system of equations is solved iteratively. 

One step of approximation is as follows: if it is 
assumed that material 2 is strongly ferromagnetic 
then Bix and Biy can be neglected with respect 
to B2X and i?2y respectively. The energy density is 
then approximated by 



/T3\ /l^z , -C ^Bx 

»(B) = +/2U'2(-^ 

^MO J2 



By „ 

h 



This function can be evaluated without solving 
an additional nonlinear equation and the magnetic 
field derives as: 
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BH-measurement for grain oriented electric steel 



Energy density for lamintated mild steel 
2 



1.5 



0.5 




Figure 4: Contour lines of the energy density func- 
tion for laminated mildsteel. 



The advantage of this approximation for FEA is 
that it can be applied without solving additional 
equations. In practical cases, however, it must be 
verified that the accuracy of the approximation is 
sufficient. 

In a first example this theory is applied to mild 
steel with an isotropic material behaviour in the 
manner of eq. [5) Figure |4] shows the lines of equal 
energy that results for a filling factor of /2 = 0.97. 
Compared with fig. [3]- the bulk material - the lines 
appear to be squeezed towards the x-axis. 



8 Grain-oriented electric steel 

Now the method of the energy function shall be ap- 
plied to model the magnetic characteristic of grain- 
oriented electric steel based on measurement. Typ- 
ical measurements are along rolling direction and 
transverse direction like in fig. [5] These measure- 
ment are performed using the Epstein-frame by 
cutting sheets into stripes along the direction to 
be measured. In these measurements the flux den- 
sity is aligned with the long edge of the stripes and 
the magnetic field strength component parallel to 
that edge is measured. The stripes are placed in x- 
direction, then the measurements gives as a func- 
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Figure 5: Measurement data for grain-oriented elec- 
tric steel in rolling and in transverse direction com- 
pared with a measurement of mild steel. 



tion B — 7> iJ that is a part of the complete three- 
dimensional function: 



H„ 



In the same way a second measurement along y 
gives an additional view. From these two measure- 
ment an interpolation scheme is established for a 
two dimensional BH-characteristic. The functions 
from measurements will be called Bq and Bqq for 
the rolling direction and the transverse direction, 
respectively. 

First the functions are integrated: 

wo{B) = / Ho{B')dB' 

WQoiB) = / HQo{B')dB'. 
Jo 



These functions are the basis for an interpolation 
rule for the energy density function w : B^ , By — > 
w{Bx, By). The main point is that the lines of equal 
energy are constructed for all values of w which 
sufficiently defines the function. For a given value w 
the inverse functions of wq^Wqq are used to find 
Bq , BgQ such that 

W ^ Wq{Bo) ^ WgoiBgo). 

The iso-line is defined to be an ellipse: 
Bx \^ f B,i 



Bo{w) 



y— I = 1 

Bgoiw)^ 



(7) 
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Iso-line for energy density function 
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Extrapolation of BH-measurement 
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Figure 6: The iso-line for energy density w is con- 
structed of be an ellipse that meats measured re- 
sults at the axes. 



Figure 7: Measurement data for grain-oriented elec- 
tric steel in rolling and in transverse direction and 
extrapolated data to capture saturation. 



Apparently the enclosed area is convex which is 
essential because a convex interpolation function 
must be found. In a practical implementation the 
flux density is given and the energy density must be 
found. A system of two equations has to be solved 
in this case. The two equations are the ellipse equa- 
tion and the requirement that the half axes of the 
ellipse correspond to the same energy. To be solved 
for Bq, Bqq: 



+ 



2 / R N 2 

B90 , 



Bx 
B~Q 

Wo{Bo) = Wgo{Bgo) 



1 



(8) 
(9) 



From the ellipse equation the equation determin- 
ing the magnetic field strength are derived by tak- 
ing the partial derivatives with respect to the flux 
density components: 



Hx 
Hi, 



1 



_B|_ 



Bx/Bq 

By/Blfi 



(10) 



Here only the two-dimensional case is shown, the 
extension to three dimensions is obvious and leads 
to a system of three equations. 

9 Application example 

The interpolation scheme shall be applied to the 
measured BH-curves shown in fig. [5] The measure- 
ments are limited to below 2 T, but when ap- 
plications in the saturation regime are targeted, 



the measured data must be extrapolated. For the 
case of the rolling direction saturation is practically 
reached and the curve can be extended using the 
permeability of one: 



Bhighfield — B sat + 



H 
Mo' 



The case for the transverse direction is more com- 
plicated. Reference [Hj) states that in any direction 
the saturation flux density is the same. Therefore 
a continuation is constructed that approaches the 
same saturation equation. This is depicted in fig. [Tj 

The application of the interpolation scheme in 
the sheet plane give the iso-line plot shown in fig. |8] 
The ellipses are clearly anisotropic but tend to be- 
come closer to circles at higher flux densities of 
«3 T when saturation is reached and the material 
becomes isotropic. 

The next step is to combine the result of the 
grain-oriented material with the lamination effects. 
Here measurements for the BH-characteristic in di- 
rection perpendicular to the sheet are needed in 
order to apply the interpolation scheme in three di- 
mensions. There are no such measurements publicly 
available. The interpolation rule with two extreme 
cases is demonstrated. One uses the BH-curve of 
the rolling direction also for the out-of-plane di- 
rection. The other case uses the measurement in 
transverse direction as BH-curve in the perpendicu- 
lar direction. The two results for the energy density 
are both displayed in fig.[9j These iso-lines show are 
cut in the plane of rolling direction and stacking di- 
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Isolines of energy for grain-oriented steel 



Isolines of energy for grain-oriented steel 
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Figure 8: Iso-lines of magnetic energy for grain- 
oriented electric steel based on measurements along 
rolling direction(a;) and transverse direction(2/); 
other points according presented interpolation 
scheme. 




0.5 1 1.5 2 2.5 3 

Figure 9: Iso-lines if energy density of grain oriented 
electric steel with lamination effect. Model 1 (red 
line) uses the BH-curve of the transverse direction 
in perpendicular direction; model 2 (blue line) uses 
the measurement from rolling direction. 



rection of the three-dimensional picture. It can be 
seen that the difference between both cases is small. 
For points with dominating rolling component the 
two models naturally give same results. With dom- 
inating stacking component the filling factor is de- 
termining the result, hence both models give the 
similar results. 

10 Conclusion 

This article proposes the use of the energy density 
function as the appropriate quantity to describe 
magnetic constitutive relations. Energy consump- 
tion is guaranteed. The dependence of the field 
strength on the flux density is found by derivation. 
The differential reluctivity is automatically sym- 
metric. 

For a thermodynamically stable material the en- 
ergy density is a convex function of the flux density. 
This leads to a positive definite differential reluc- 
tivity and the finite-element-method is stable. 

The method is in particular suitable for aniso- 
tropic materials. The BH-law for anisotropy from 



lamination is demonstrated. For grain oriented steel 
an anistropic material law is derived from two mea- 
surements (along rolling direction and transverse 
directions). The interpolation for arbitrary direc- 
tions is done using the energy density and by build- 
ing isosurfaces for this function by using measured 
date as starting points. This BH-law is used as in- 
put to the previously found lamination law. 

If measurements are available in intermediate di- 
rections the data can included in the isosurface con- 
struction. 

It must be repeated that be presented method is 
only applicable for anhysteretic material because 
only then the energy density is a function of one 
variable (the flux density) alone. Furthermore 
the measurements must be performed such that 
the material is always in a thermodynamic stable 
state, on then the convexity condition is valid. 
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